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Методи анал1зу сигналав звукового д1апазону 


масивом сенсор1в моблльного робота 


В статье рассматриваются теоретические и практические вопросы построения системы пеленгации (т.е. 
определения направления на источник в трех измерениях) слабых звуковых колебаний в воздушной среде. 
Результат получают, основываясь на усредненных во времени полиномах второй степени относительно 
акустического давления и/или компонент колебательной скорости. 
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У стати розглядаються теоретичн! 1 практичн! питання побудови системи пеленгаци (тобто визначення 
напряму на джерело в трьох вимтфах) слабких звукових коливань в повтряному середовици. Результат 
отримують, трунтуючись на усереднених в час! полномах друго! мри вдносно акустичного тиску Габо 
компонент коливально! швидкост!. 
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Введение 


Пусть на фоне шума окружающей среды имеется слабый (например, на 10 дБ 
меньше в точке приема) сигнал сосредоточенного малошумного объекта. Комбиниро- 
ванный векторный приемник измеряет в полосе частот анализа акустическое давление Р 


и колебательную скорость Г. Обрабатывая измеренные величины, требуется обнару- 
жить наличие объекта и найти его пеленг. Известно решение, основанное на использо- 
вании вектора потока мощности и спектральном анализе. Однако эти решения трудно 
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применимы для реализации на микроконтроллерах, находящихся на борту мобильного 
робота, выполняющего поиск источников сигналов. Ниже предлагается более реализу- 
емый и не менее эффективный подход к решению такой задачи. 


Основные соотношения 


Введем обозначения. Плотность среды р, скорость звука с, компоненты колеба- 
тельной скорости в декартовой системе координат И,Г,,Г,, оси 1 и 2 направлены 
горизонтально, ось 3 направлена вертикально вверх. Направления осей | и 2 считаются 
известными, например, определяются с попощью компаса или СР$ системы. Обозна- 
чим И, =Р/Хрс) . В полосе частот шириной Р [Гц], полагая шум окружающей среды и 
шум объекта близкими к белому, независимые измерения сигналов проводятся с 
периодом т=1/(2К)с в моменты времени &,1= (1,2,...М). Угловыми скобками будут 
обозначаться усреднение во времени, например: 


са [иезиль м 


Для величин <У\/. >,< И, > ‚ подразумевается то же. Диаграмма направленности 


излучающего шум участка поверхности [1, с. 659] ©*(©) = Г, соз”"(О), где Г, — интен- 


сивность излучения при © =0,а © - угол отклонения от вертикального направления. 
Известно, что интенсивность принимаемого сигнала от участка поверхности пропор- 
циональна его площади и обратно пропорциональна квадрату расстояния. Число т 
лежит [1, с. 663] в пределах от 0,5 до 1 для частотного диапазона до 300 — 400 Гц и 
т=1, 2, 3 для диапазона частот выше 500 Гц, с увеличением скорости ветра и с 
увеличением частоты число ти растет. В [2] показано, что число т может быть равным 
0,5 даже до частоты 700 Гц. Можно обозначить 


СЕТ. 

Путем несложных вычислений можно получить: 
<И, >=-Т@2т+2)/(2т+1), 
<И; >=Т,(т+1)/т=Е,, 
рай ЕТ: 
Тогда для целей моделирования сигналов можно использовать четыре независимых 
случайных процесса с нулевым средним и единичной дисперсией &,(1,) и 5,(1,): 

а ЕН 

Шум окружающей среды можно смоделировать следующей системой: 


т 
ЕЕ т (1) 
2т 1 
У =-=.. Е +5. .|Е : 2 
` аж, | Эт (2 


1 

и = 1 | ее (3) 
| 1 

Г, ==, |Е; 2(т-+1' (4) 
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Далее необходимо ввести систему ограничений вида: 
2(т+1) 
2т+1' 


Г, =И +: › где 4= (5) 


И =< ИИ > ре, 

И, =< ИУ, > ре, 

И! =<ИИ > рс, 

И’, =< И, > рс. 
Для объекта, излучающего шум в горизонтальном направлении, вектор потока 
мощности с компонентами И, И, и вектор с компонентами Й,, И’, при болыпом 


времени усреднения совпадут. При конечном времени усреднения флуктуация век- 
тора будет существенно меньше, т. е. эффективное отношение сигнал/помеха станет 
лучше (примерно на 3 дБ, при т =1,5). Сигнал малошумного объекта можно описать 
выражениями вида: 


Р=0, 


Г, ==. \/Ё., 
И =-=.\/Е, с0$(ф), 
Г, =-=.\/Е. зи (ф), 

где фФ — угол пеленга объекта в горизонтальной плоскости, #,(1) случайный 
процесс с нулевым средним и единичной дисперсией, а величина Ё, характеризует 
интенсивность излучения шума объектом. Обозначим: 

ЙЕ, 
О, =2<ИГ, >. 

Из этих уравнений по полученным в результате измерений величинам И и И, 
можно вычислить величины Ё, и Ф (с точностью до слагаемого, кратного п ) и можно 
построить вектор И’ с компонентами 

| = Е. со5(ф)рс, 
И, =Е, зш(ф)рс. 

Направление вектора Ц’ определено с точностью до смены знака. Смысл вектора 
такой же, как и вектора потока мощности У или вектора \" с уменьшенной флуктуа- 
цией. Шум окружающей среды искажает измеренные значения векторов И и (', но 
согласно формулам (1-4) шумовая составляющая, вносимая в вектора И” и И’ прак- 
тически независима друг от друга. В вектор Й” дают вклад произведения < &,Е, > и 
<=)Е, >, аввектор И’ дают вклад произведения <&5, >, <5^>, <,” > и, следова- 


тельно, взвешенная сумма этих векторов будет зашумлена меньше, эффективное от- 
ношение сигнал/помеха возрастет. Другими словами, вычисляются два вектора: 


М, = (а, + ВИ) (а + В), 
М, = (И, + ВИ, (а + В), 
М, = (а, - ВИ) (а + В), 
М, = (а, - ВЦ.) (а В). 
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М - с компонентами М.и М, иМ ' - с компонентами М, и М,, и в качестве 


ответа выбирается больший по модулю. Таким образом, вектор потока мощности по- 
зволяет нам правильно выбрать знак вектора (У, а вектор С” вместе с вектором И” 
позволяют снизить флуктуацию результата. 


Выводы 


Коэффициенты @ и В могут быть выбраны в зависимости от желаемой цели. Таких 
существенно разных целей видится две, а именно: точно найти направление на объект и 
точно определить уровень его сигнала. Для мобильного робота, работающего в режиме 
поиска, актуальными являются обе задачи, но это не означает, что их надо решать одно- 
временно. Пусть на данном этапе актуальной является задача наиболее точно 
определить направление на объект — источник сигнала. Тогда в соответствии со 
структурой системой [3] для каждой пары сенсоров (микрофонов) находятся вектора 
М‚и М, иМ' - с компонентами М! и М, , из пар выбирается больший по модулю и 


эти вектора собственно и будут определять конус направления на источник сигнала. 
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